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Abstract. Turbulent hydrodynamics have been in the focus of interest of statistical physics at least since 
69 years, when Kolmogorov's famous 1941 publication set the basics for a framework of research. The 
Renormalization Group is since a long time believed to be a very promising tool for the analysis of 
turbulent systems, but previous attempts had relatively little success. In this work, we present an overview 
over the problems involved, as well as our approaches to overcome them, and first, promising numerical 
results. 

PACS. 47.27.ef Turbulent flows - Field-theoretic formulations and renormalization 



1 Motivation 

Since the first work on statistical hydrodynamics by Kol- 
mogorov [20J, turbulence basically remained an unsolved 
problem of classical physics, even though the fundamen- 
tals seem to be simple enough - the Navier-Stokes equa- 
tions just express the conservation of momentum of an 
incompressible fluid. 

■ Nevertheless, experimental results on the moments of ve- 
' locity differences (the structure functions to be defined 
, below) have yet to be understood on a theoretical basis. 
• While Kolmogorov's assumptions of scale-independence 

fail, the intrinsic scale-dependence that is responsible for 

■ the intermittent exponents have not been deduced from 
' first principles so far. 

, A promising approach seems to be the Renormalization 

■ Group (RG), which aims to describe the dependence of 
I the correlation functions of a given field theory on the 

■ scale on which the system is observed. Beginning with the 
' works of Forster, Nelson and Stephen [TT], numerous at- 
, tempts have been made to apply the various formulations 
' of the RG to turbulent hydrodynamics, but until today 
] the observables proposed by Kolmogorov could not be de- 

■ duced in accordance with experiment. For some work on 

■ the RG approach to turbulence, see [51II] . 

The formulation of the RG most suitable for the study of 
turbulence is, in our opinion, the so-called Exact Renor- 
malization Group (ERG), due to Wilson [5^1155] . Although 
the distinction between "field theoretic" and "exact" RG is 
merely artificial, the latter explicitly involves the generat- 
ing functional of the theory to be studied, which can be 
simplified and re-formulated to suit the analytic methods 
involved. It is especially helpful, as we shall see, for the 
analysis of a theory with constraints, like in our case the 
incompressibility condition. We will show how to apply the 
well-known Faddeev-Popov gauge-fixing method to work 



out the functional integral measure. In this step we differ 
from previous work, which in one way or another omitted 
the incompressibility condition and/or the pressure term 
of the Navier-Stokes equations. 

The RG-flow, as we shall discuss in detail, can be under- 
stood as the continuous way of calculating all Feynman 
graphs of the theory. Keeping this in mind, we established 
a numerical algorithm that simulates the RG-fiow by inte- 
grating out the corresponding graphs. We therefore take 
advantage of the freedom of choice of a cutoff-function for 
the propagator. We end up with lengthy, but straightfor- 
ward rate equations, that can be iterated quickly and up 
to a high number of involved field operators. 
We will also show that the predictions of Kolmogorov 
can be found as the trivial scaling solutions of this the- 
ory, and that non-trivial structures in coupling space ex- 
ist. We have tested the algorithm on well-known theories; 
the identification of intermittent exponents in turbulence, 
however, has not yet been accomplished. 

2 Turbulent Preliminaries 

We shall introduce the basics that are needed in this work, 
and sum up Kolmogorov's predictions from 1941 (K41) as 
far as we aim to falsify them, rather than to give a full ac- 
count of Kolmogorov's theory. Reviews can be found for 
example in [^[Bpgpi]- 

2.1 Navier-Stokes Equations 

We shall mainly work with the full Navier-Stokes Equa- 
tions (NSE) given by 

dtVa + VpdpVa - I'V^Va + -daP = 0, (1) 
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wliere v denotes tlie _D-dimensional velocity field (in 
Navier-Stokes turbulence, D is either 2 or 3), v the kine- 
matic viscosity, p the scalar pressure field and p the density 
of the fluid. Eq. ([T]) expresses the conservation of momen- 
tum in an infinitesimal volume element of the fluid. As 
these are three equations for four degrees of freedom, it is 
clear that we need an additional constraint, determining 
the pressure field. In incompressible turbulence, this con- 
straint simply states that the velocity field is divergence- 
free, 

M{v) dcVa = 0. (2) 

Starting from any given initial condition, the velocity field 
obeying Eq. ([T|) has to develop into a constant field, as the 
dissipative term converts more and more energy into heat. 
If we aim to understand fully developed turbulence, sta- 
tistically homogeneous in space and time, and statistically 
isotropic in space, we need a mechanism to insert energy 
into the system, so that an equilibrium flow can develop. 
The standard way of providing this is to add a stochas- 
tic force (stirring force) fa to Eq. ([1]) that is long-range 
correlated: 



notes the fundamental solution of the Laplacian, e.g. in 3 
dimensions: 



1 



Awlxi - X2\ 



(6) 



Different forms have been tried for F, though in the con- 
text of the NSE it is widely believed that the form of the 
stochastic force does not influence the statistical charac- 
teristics of turbulence. On the other hand, it should not 
be forgotten that e.g. in Burgers turbulence, the intermit- 
tent exponents (to be deflned below) clearly depend on 
the choice of the stirring. 

Eq. (2)) is sufficient to eliminate the pressure term, as can 
be seen when deriving the solenoidal NSE, as we shall do 
in the following. In a first step, let us operate onto Eq. ^ 
with a divergence operator. Then the first and the third 
terms drop out, as the field is divergence free: 



dtVa + Vpdf^Va 



-daP = fa- 
P 



(3) 



We then invert the Laplacian, arriving at 



The idea is to bring energy into the flow on large scales, let 
large structures decay freely into smaller ones until the en- 
ergy is finally dissipated into heat (Richardson-cascade). 
We model the stochastic force to be Gaussian distributed, 
with (5-correlation in time, and a long-range correlation 
function in space: 



1 

-P 
P 



dp 
V 



dp 

V 



(7) 



(8) 



{f{x^,t)f{x2,t')) 



F-\xi,t,X2,t) (4) 
-eSit-t')V-^ixi,X2), (5) 



where e is the local energy dissipation rate (not to be 
confused with the parameter of the e-expansion). V^^ de- 



which is the above mentioned condition for the pressure 
field. 

One might ask whether the inversion of the Laplacian 
leads to a unique solution for p. Now, two solutions might 
at best differ by a harmonic function, which is either con- 
stant or growing without limits. The second option is not 
physical, the first one not relevant as we are only working 
with pressure differences. 

To obtain the solenoidal NSE, replace the "solved" pres- 
sure field into Eq. ([3]): 



dtVa - vV'^Va + ^a 



dadf, 



[v^d^Vp) = Sap 



dad, 



V2 



dtVa - ly'^^Va + Papiv-fd^Vp) = Papfp- 



(9) 

(10) 



Observe that the operator Pap ~ Sap — is precisely 
the transverse projector known from electrodynamics; a 
feature that we will use later on. 

From now on we shall investigate the solenoidal NSE. It is 
important to keep in mind that these are only equivalent 
to the full NSE as long as incompressibility is ensured. 
Also observe that Eq. (jTU]) is non-local, as it involves the 
inverse of the Laplacian operator, the integral kernel of 
which is of the form (|6]). 



2.1.1 Gauge-lnvariance 

As a preparation for a Faddeev-Popov method, that we 
will need in the following section, we shall now rewrite 
the formula in a gauge-invariant way. As the transverse 
operator P projects a field onto its incompressible parts, 



and knowing that the fields we are interested in are trans- 
verse a priori, it is easy to see that 



PapVp = Sa 



dadp 



Vp = SapVp = Vo 



(11) 



SO that we are free to replace v by Pv in Eq. ([3]): 

N(i;) := dtPapVp - vV'^PapVp + PapiPjSVsd^PpeVe) 



Papfp- 



(12) 



Now, even if the resulting equation looks more compli- 
cated, it is invariant under the same gauge transformation 
known in electrodynamics. 



Va ^ Va + daA{x) —: U{v). 



(13) 



Constraint ^ is still required, but it now acts as a gauge 
fixing term. 
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2.2 Structure Functions and Intermittent Exponents 2.3 Burgers Equation 



In 1941, Kolmogorov published his famous article, intro- 
ducing a statistical framework for turbulent hydrodynam- 
ics |20J. As the theory is Galilean invariant, he proposed 
that the observables should be functions of velocity dif- 
ferences, and more specific, he considered the so-called 
velocity increment 



;(r; x) = {v{r + x) — v{r)) 



(14) 



the difference of the velocities at two poiirts separated by 
a vector a;, projected onto the unit vector in x-direction. 
Suitable observables appear to be the structure functions 
of order p. They are defined as the p-th moment of the dis- 
tribution of the absolute value of the velocity increment: 



Sp{x) {\vinc{r;x)\'')r 



(15) 



The average is taken over all points r of a realization of 
the turbulent flow. In the case of homogeneous turbulence, 
we suppose this to be equivalent to an average over all 
histories v{x,t), which we will perform in the functional 
integral. 

Kolmogorov proposed the existence of a smallest leirgth 
scale A, also called the "dissipation scale", below which 
physics is no longer dominated by turbulence, but by dis- 
sipation. Dimensional analysis leads to 



(16) 



where e is the (constant) dissipation rate. Assuming that 
the turbulent cascade of decaying vortices happens on 
scales much larger than A, we might hope that observ- 
ables doir't depeird on it and are thus self-similar, which 
means power-law functions of the scale: 

Sp{x) oi {ex)^" . (17) 

Again, Kolmogorov deduced by dimensional analysis that 



P 



(18) 



It has long been stressed that the fundamental assump- 
tion, namely the indepeirdence of the smallest scale A, is by 
no means natural, and is in general not fulfilled in critical 
systems. This could lead to a scale dependent dissipation 
rate (or viscosity) and the breakdown of scaliirg law (fT7|) . 
Even though general agreement on this point seems to be 
common, the scale dependence could not yet be deduced. 
In case that a typical (macroscopic) length scale L can be 
identified in the system, the Reynolds number is defined 
as 

\ 10/3 

A 



Re 



(19) 



L might be the radius of an obstacle of the flow, or, more 
useful for our purposes, the correlation length of a choice 
of the stochastic force. Eq. (|19p coincides with the more 
common definition 

Re = — (20) 



if the typical velocity U is defined to be 
U = (eL)i 



(21) 



As we shall refer to it later, we mention Burgers' equation 







(22) 



as a simpler model for turbulence, which can be seen as 
fully compressible hydrodynamics. Further constraints are 
not needed as a pressure term is missing. Burgers turbu- 
lence, or Burgulence, is investigated in one to three space 
dimensions (see e.g. |t3j). It is local, and its solutions typ- 
ically contain structures that resemble the fundamental 
solution of the Hopf-eq. 



dtVa + VfidpVa = 0, 



(23) 



as shown in Fig. ([T]). These are regularized by the viscos- 
ity term to a finite step width (similar to the behaviour of 
the KPZ-equation, see [21]). 



Shock structure in Burgers' configuration -|- 



-I- 
+ 



+++" 



+, ,++ 



Fig. 1: Typical solution of Burgers' equation from a 
Monte-Carlo simulation [lOj. The graph shows v{x, t) as a 
function of x in periodic boundary conditions at constant 
t. 



As published in a different article [lOJ, we investigated 
a generating functional of Burgers equation by means of 
Monte-Carlo methods. These investigations have shown 
that the generating functional formalism is suitable for 
the study of turbulent systems, even if almost-singular 
structures are involved. 



2.3.1 Exponents for Burgers Turbulence 

In general, the intermittent exponents of burgulence 
strongly depend on the form of the stirring force; for an 
overview see e.g. [14]. In the canonical case, that of freely 
decaying burgulence in 1+1 dimensions, analytic calcula- 
tions lead to 

= min {p, 1) . (24) 

This scaliirg can be interpreted as the signature of a bifrac- 
tal decay channel. 
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3 Generating Functional 

The basis of any analysis using the Exact Renormalization 
Group (ERG), first introduced by Wilson, or many other 
field theoretical methods, is the geirerating functional. For 
turbulence several attempts have been made to define it 
as a functional integral, all in one way or another omitting 
the condition of incompressibility. We shall begin with the 
derivation of the Martin-Siggia-Rose functional [23] for the 
solenoidal NSE, and than show how to respect condition 
([2]) by means of a Faddeev-Popov procedure. 



3.1 Fine-Grained Distribution 

Our starting poiirt is the so-called fine-grained probabi- 
lity distribution for the field u, obtained by counting all 
possible solutions to the NSE. 



The (5-function is now written by means of a functional 
Fourier-transformation, introducing a new formal (book- 
keeping) field u: 



S{N{v) -Pf)(x / Vu e'(".N(^)-P/) 



(29) 



If we formally define the action S'l by 



J { dvf)(x') J 

(30) 

where the determinant still has to be evaluated, we can 
identify the elements for our Feynman-rules. For the time 
being, we want to focus the attention on two parts, which 
together lead to the famous 6'(0)-problem: 



Z[J] 



Du(^(5(u-N~i(P/))e('''"')^ , (25) - w-(diffusion)-propagator: 



where N is defined in Eq. (fT^)P 1. and we defined the ab- 
breviation 



U 



V 



d xid X2dtidt2 X 
y.^i{xi,ti)A{xi,ti;x2, t2)v{x2,t2). (26) 



A few remarks are in order: 

— Of course, N^^ is not to be understood as the (not 
necessarily existing) inverse of an operator, but as a 
multi-valued operator counting any solution v for a 
given realization of the random force /. Observe that 
we are working with a functional 5-function, meaning 
that we are searching for histories v{x, t) that solve the 
NSE for all x and rather thair a realization v{x,ti) 
at a given time ti, depending on some initial condition. 

— The average (•) is to be understood as an average over 
all realizations of /, replacing the spatial average in 
(|15p. It is the basic assumption of this work that for 
homogenous and isotropic turbulence, these averages 
are interchangeable; we are supported by our results 
for Burguleirce published in |10| . 

Making the average over all realizations of the stochastic 
force explicit yields 



Notice that the corresponding bare two-point-function, 
also called response function, is proportional to the 
Green's function of the diffusion equation, applied to 
transverse fields: 



{u{x2,t2)v{xi,ti)) OC 



1 



(31) 



We have to choose the retarded Green's function to 
ensure causality of the theory, so 



{u{x2,t2)v{xi,ti)) OC 6{t2 - ti). 



(32) 



uvv-veitex: 



Z (X J J 2?/(5(u-N"i(P/))e" 



|(/,F/)+(«,J) 



(27) 



U 



Now the argument of the (5-function is multiplied by N, 
which leads to a functional determinant that we shall dis- 
cuss in the next paragraph: 



^ V 



N V 



hU-Ff) 



I 5vfi{x') 



(28) 



In a graphical expansion, it is required to calculate the 
following loop: 



^ At this point, we adopt the notation of L'vov and Procaccia |22| . 
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V 




where 



-z(u,NH-P/)--(/,VV) 
+*(^*'^^)- (39) 



Fig. 2: u-w-loop of the interaction term 



As the retarded Green's function of (not only) the dif- 
fusion equation is proportional to 9{t2 — ii), this loop 
is proportional to the seemingly ambiguous quantity 
9{Q). We will show in a following paragraph how this 
is solved by a choice of discretization of the functional 
integral. 

To illustrate what the above action actually implies, we 
integrate out the non-physical fields / and u for a moment. 
Observing that first the integration over /, and then the 
one over u are Gaussian, we can simplify Z to: 



Z 



with 



Vv e^^^^I^ldet 



1 



\ 5n^{v){x) 



S2M = -(NH,FNH). 



(33) 



(34) 



We see that field configurations not solving the NSE are 
admitted, but suppressed by a Gaussian weight that we 
can control numerically, and that can be related to a phys- 
ical constant by Eq. dS]). 

This expression is already a generating functional, as we 
are able to extract all correlation functions using func- 
tional derivatives. Before we can start to simulate RG- 
flows, we still have to bring the determinant into a suitable 
form. 



3.2 Functional Determinant 

A straightforward way of writing the functional determi- 
nant is by using ghost fields. 



det 



6vp{y) 



(35) 



Vrv^ g-*(^-,(9*p-.v^))^+»(>^Mi>^)^ (36) 



where we dropped a field-independent term and defined / 
to be the non-linear part of N, 



Ia[v] = VpdfiVa- 

We thus end up with the functional 



Z = VvVuVf V^* e 



(37) 



(38) 



It is important to keep in mind that the ghost fields, 
though anticommuting with each other, don't pose any 
particular problems in the numerical procedure. In prin- 
ciple, it is very much feasible to take them into account, 
and we performed several RG-runs while keeping track of 
the ghost terms. 

The determinant has a simple graphical interpretation, as 
we see directly that it exactly cancels out the M-i;-loop 
shown in Fig. ©: From eqs. (PT|) and we see that 

the and the uw-propagator are identical. The reader 
can easily check that ip* , leads to two terms similar 
to uN[i;], but with u replaced by "0* and one w-fields re- 
placed by -if). When this vertex is closed to a loop by means 
of a V'*'0-pi'opagator, this is numerically identical to graph 

This greatly simplifies matters for us, as our numerical 
program sorts and calculates contributions to the RG- 
flow according to their graphical representation. Rather 
than simulating two additional fields, and calculating all 
the graphs, we are thus allowed to drop a certain class of 
graphs. The cancellation of certain averages can be proven 
even non-perturbatively, using the BRS-invariance of the 
action. 



3.2.1 BRS-lnvariance 



If the functional determinant is expressed in terms of 
ghost fields, this yields an extra symmetry, also called 
BRS-invariance. We observe that the action ([M)) is indeed 
invariant under the following infinitesimal simultaneous 
changes: 



5u = 0, 
5il)* = 0, 
5v = eV*: 
(5-0 = ieu, 



(40) 
(41) 
(42) 
(43) 



that for obvious reason have also been called "half a super- 
symmetry". From the Ward identities of this symmetry, 
the desired result follows on a non-perturbative level: 



ov 



(uN) . 



(44) 



Both sides of this equation can be interpreted as the sum 
of the corresponding graphs as explained in the previous 
section. For details we refer the reader to the explict proof 
in 1271. 
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3.2.2 The 6I(0)-Problem 

We have seen in (|3.ip the typical appearance of the seem- 
ingly arbitrary number 6{0) in the retarded Green's func- 
tion of the diffusion equation. It is of course linked to 
the Ito-Stratanovich-dilemma; Zinn- Justin [35J shows ex- 
plicitly how it is connected to the problem of operator 
ordering. To clear the mist, we want to show in this sec- 
tion how 6(0) is fixed by the choice of discretization of the 
functional integral, namely the time derivative. 
In any case the 0-function is required to obey 



(x) =0, Vx<0, 
(x) 1, V a; > 0. 



(45) 
(46) 



We further require the fundamental theorem of calculus 
to be valid, so 



e{o) 



dt9{t) dt. 



(47) 



We now have to define a differentiation. We distinguish 
different choices by a parameter a, suggesting 



dtf{t)^yiwJ-[f{t 



1 



lit 



1 



r) , (48) 



with — 1 < a < 1. Notice that a = corresponds to the 
symmetric, and a = — 1 to the pure backward derivative. 
We conclude 



^(0) = lim i f \e{t 



a + l, a— 1,., 

r) - 6{t + ^— r) \ dt. 



(49) 

The term involving the second 0-function does not con- 
tribute since t < in the integral. The first ^-function 
only gives a contribution for 



t > — r. 



(50) 



Incompressibility doesn't matter. Though not stated 
explicitly, this appears to be the point of view taken in 
those works that neglect the additional condition ([2]) 
without any further comment. This is more than ques- 
tionable, because there is a whole class of equations 
that only differ by the compressibility condition, and 
give different statistics. An obvious example is Burgers' 
equation (|22p . which models fully compressible fiuids 
and shows bifractal scaling of the structure functions. 
Incompressibility is enforced by a choice of bound- 
ary conditions. This argument goes as follows: In the 
solenoidal form, it can be clearly seen that an incom- 
pressible fiow stays incompressible even without en- 
forcing it by a particular equation, as the random force 
term on the one hand, and the former pressure term 
on the other hand lead to incompressible contributions 
to the fiow. On the other hand, all compressible parts 
of a given flow tend to die out due to the dissipation 
term. This leads to the observation that condition ([5]) 
is always ensured if we impose an incompressible fiow 
for t — >■ — oo. This observation is indeed useful for some 
calculations, e.g. direct numerical simulations, as they 
have a definite boundary at t = 0. In our case, we aim 
for statistics over whole histories of the fields, for all 
times - and in negative t-direction, as the argument 
itself states, any compressible part of the flow is going 
to be amplified. Now the slightly compressible fiows 
and the incompressible ones lie dense to each other 
in functional space, so that in any numerical applica- 
tion we would lose control of the boundary conditions 
completely. 

Incompressibility is ensured by some standard integra- 
tion measure. What is meant here is that incompress- 
ibility is considered to be included in an integration 
measure "Dvinc (we are going to proceed in the same 
way in the next paragraph), but then the measure is 
treated as a standard Lebesgue-like measure. This is 
certainly not legitimate; e.g. an integral like 



3-(finc,At)inc) 



(52) 



We thus find 



1 

61(0) = lim - / dt = 



a + 1 



(51) 



Thus 9{0) is fixed by the choice of the time derivative. 
Notice that for the symmetric (Stratanovich) derivative, 
a = 0, we find 6'(0) = i, while in the pure backward (Ito) 
case we get 6{0) = 0, as expected. 



looks like a Gaussian integral, but will in general not 
be so. 

We thus conclude that incompressibility should well be 
taken care of. It might be surprising that the way to do 
this is rather straightforward and well known for gauge 
invariant field theories. 



3.3.1 ,, Strong" attempt using (5-function 



3.3 Incompressibility Condition 

Before we demonstrate how to derive a generating func- 
tional respecting condition ^ for the full Navier-Stokes 
equations, let us briefly discuss three common positions 
concerning incompressibility: 



We implement condition ^ through a second functional 
Fourier transformation of a 5-function, 



(53) 



The advantage of this is that the generating functional is 
exact, apart from the definitions of some constants. The 
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major drawback is that we will have to apply the Renor- 
malization Group in a discretized way, which requires an 
approximation of the generating functional. Now a con- 
straint like ([551) can, after approximation, lead to a field 
theory that has no solutions at all, as it is driven from the 
manifold of solutions further away by every iteration of 
the RG-process. 

From a technical point of view, the additional field makes 
simulation times longer and most calculations more elab- 
orate, but does not represent qualitative problem. 

3.3.2 Faddeev-Popov Method 

A proper way to treat the compressibility condition is by 
means of the Faddeev-Popov method. Let us first general- 
ize our considerations to the case daVa = c. We multiply 
the functional by 1: 

1^ JvUS{U{v)-M-\c)) (54) 

= /pL/^(M(t/(.))-c)det{M^}, (55) 

where U denotes the gauge transformation as defined in 
(fT5|) . The determinant can be written explicitly as 

dct{«^}^det{«M| (56) 

^det{^iM^^|±M))| (57) 

= dct{v2}, (58) 

and can be dropped as it is independent of any field. 
We now multiply the functional by another (field- 



We discarded convention ([2^ because in the next section 
we shall directly manipulate some of the terms that would 
otherwise be hidden by the notation. 



3.4 Galilean Invariance 

A remark concerning Galilean invariance, as analyzed by 
Hochberg and Berera[5], is here in order. The path inte- 



independent) factor, namely 

e-^'-'^^'^l (59) 

Now our original action is gauge- invariant, meaning that 
the generating functional does not depend on the value of 
c. We might choose a particular one, or as well integrate 
over all possible c: 

Z oc J VcDuVvVf e"^'t''''''^'<5(M(w) - c) x 

xdet|^^^44^1e-^(-). (60) 

The reason to integrate over all c is to evaluate the 5- 
function, finally arriving at 

Zoc|p.P.P/e-^["-/Uet{^^^^m| (61) 
with 

S[u, V, /] = ^{dv, dv) ~ i{u, NH - Pf) 

-^(/,VV)- (62) 

The gauge fixing term appears as a Gaussian distribution, 
and only in the limit k — > incompressibility is enforced 
strictly in the integrand. This method overcomes the dis- 
advantages of the (5-function method, as deviations from 
solutions aren't completely forbidden, but only suppressed 
by a Gaussian weight. The limit can be taken when re- 
sults have been obtained. We shall from now on work with 
the functional (j62p . but formulated in wavenumber space, 
which is: 



(63) 

gral, like the NSE, is invariant under the changes 

v{x,t) ^v'{x\t')+c, (64) 
X = x' + ct, (65) 
t = t'. (66) 

In any attempt to calculate averages or correlation func- 
tions of the field itself, this might be a problem that could 
be overcome by another application of the Faddeev- Popov- 
method. In our case, however, it is of no practical con- 



S[u,V,f]=J ^--jS{p + q)dt 

~iUa{p)vq'^Pal3Vli{q) +iUa{p)Pafs{q)fl3{q) + ^fa{p)q'^ fa{q) 

Ua {p)PaB (p) {P-,5 {q)vs {q)r^ Pae {r)v^ (r) ) . 
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sequence, as we are aiming for the averages of velocity 
differences, which are gauge invariant. 



3.5 Non-local Interactions 

The action we derived still contains interactions of the 
type 

UaPap{PfSVsd^PpeVe), (67) 

which are non-local in physical space, but of a very sim- 
ple form in wavenumber space. Independent of the space 
where the theory is formulated, we need to rewrite (|67p iir 
a local way: In physical space, non-local interactions are 
at best cumbersome; in wavenumber space, we are going 
to sort the terms of the action according to their power of 
momenta, so we try to avoid ^--interactions. 
We shall proceed in two steps: we will first re-define non- 
physical fields, and then introduce new fields to transport 
the non-locality of interactions. We will end up with a 
lengthy, but local action that suits our needs for further 
analysis. 

To make the non-local irature of the iirteractions more ob- 
vious, the functions within this paragraph are functions 
on physical space. 



3.5.1 Field Redefinition 

If we redefine the non-physical fields in the following way: 



and introduce 

we can rewrite the action as 



(68) 
(69) 
(70) 



(71) 



iUaPapiP'fSVsd-^Pf^eVe) ^ -iQapUad^Pfi^V^P^sVS, 



(72) 



(73) 



(74) 

The functional determinant of these transformations is 
field independent and can thus be omitted. 



3.5.2 Introducing new fields 

Some non-local terms of the general form 



(75) 



for some fields K and L still survive. We are able to solve 
this problem by introducing new fields. These can be un- 
derstood as transportiirg fields that "carry" the non-local 
interaction from one place to another, thus replacing it 
by two local interactions and a propagator. Formally this 
means inserting a Gaussian integral of the type: 



VMe 



(76) 



where we define M to be a new field. This leads to a new 
kinetic term i(M, V^Af) in the action, which is indepeir- 
dent of all physical fields. We now shift the variables as 



1 



1 1 



M := M + X-'—K + \ -L. 



2 V2 



(77) 



The poiirt is that 



- MV'^M + K—L = -MV'^M ~ 2\-^MK - XML 
-X-'K^K Ix'L^L. (78) 

So we get rid of the original, noir-local interaction (|75p by 
replacing it by a new kinetic term for M and new inter- 
actions. Two of them are still non-local, but of the much 
simpler type 

1 



yL—L, 



(79) 



for example. We can deal with them by the same method 
to get a local action finally. We add again a Gaussian in- 
tegral for a new field, say N, and define 



from which we end up with 



1 



1 



L = -NV^N - iXNL. 



(80) 



(81) 



The constant A is needed so that the new fields get a def- 
inite dimension. 

This procedure might appear unusual, but is in fact quite 
commoir. The reader may think of electrodynamics: Wheir 
a light is shining and is seen from a distance, this might 
be called a non-local interaction. In theory, a local elec- 
tromagnetic field is introduced, that propagates from the 
source to your eye. Now, the equations of motion for the 
electromagnetic field turn out to be fairly simple, so that 
we are used to think of the field as being fundameirtal. In 
the present case, think of the fields and M as merely 
helpful auxiliary quantities, though it might in other cases, 
like Burgulence, be advisable to model structures as irew 
fields explicitly, and observe their interactions in an effec- 
tive field theory. 

Applying the method discussed above to action we 
arrive at the local action 
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^loc = Oo[v,U,f] + Oi[(/)i,(/)2,03,/, U] +O2[0\(/)2,(/)3,/,'y,'«], 



(82) 



with 



On = 



Oi 





/ d^q 






d^p \ 


f d^q 


(27r)^; 





S{p + q)dt - iUa{p)dtVa(q) + iua{p)fa{q) 



S{p + q)dt[ - X(j)i{p)qafa{q) + 2X ^ 4'2{p)<laUa{q) 



+ 2iX ^(piipjqaUaiq) - i\(j)z{p)qafa{q) 



d^p \ f d'^q \ / d^r 



{2n)D J \{27:)D J V(27r) 



S{p + q + r) dt Ua{p)vi3iq)rpVa{r), 



(83) 



O2 



d^p \ ( d^q 



{2tt)D J V(27r) 



S{p + q)dt - wua{p)q^va{q) + ^^(t>k{p)q^(t>k{q) 



+ ■^faip)q^fa{q) 



d^p \ ( d^q 



{2Tr)D J \{2tt)D J \(2 
>^Pa4'3{p)vi3{q)rpVair 



fe=l 



J 5{p + q + r)dt I - i\pa4>i{p)vj3{q)ri3Va{r) 



(84) 



Here, (f>i,(t)2,(f>3 denote the new scalar fields. The terms 
are already sorted by their order of derivatives. 
This is the result for the action S. Depending on how the 
determinant (pSj) is expressed, other non-local interactions 
may have to be rewritten in the same way. 



3.5.3 Discussion 

In this section we have shown how to transform non-local 
into local interactions: either by redefinition of unphysical 
fields, or by introduction of intermediate propagators. A 
major drawback will be that we have to approximate this 
action to account for RG-transformations, and a common 
way is the derivative expansion. Now, due to our "local- 
ization" of the action, the original terms have been mixed 
concerning the order of derivatives. Moreover, the number 
of derivatives has increased for most interactions, which 
means that we would have to expand the RG-flow to a 
high order in the derivative expansion. Also, the number 
of fields involved makes it a cumbersome work to do, even 
to the lowest orders. 

Nevertheless, this expansion is feasible to any order in 
derivatives, as shown in [18] for the first two orders. The 
results are rather lengthy rate equations, which will not 
be elaborated on here. 



4 Renormalization Group Considerations 

The renormalization group has its origins in the work of 
Gell-Mann and Low in 1954 [150. The first works concern- 
ing the Wilsonian Exact Renormalization Group (ERG) 
have been published in the early 70s [5^155] . based on 
Kadanoff's block-spin picture (TH]. It is surprising that 
some very basic questions, e.g. concerning the renormal- 
ization step and the anomalous dimension, are still being 
discussed. We will consider this point in detail in para- 
graph 14.41 as it seems to be in order (at least to us) to 
address these matters, especially the anomalous dimen- 
sion and the graphical representation of the flow. 
In this section, we remind the reader of the origins of the 
ERG, and wiU try to give a basic understanding of the flow 
equations. We will not repeat the derivation of the equa- 
tions, as this can be found in a number of articles, but 
quite a bit of the graphical representation of the differ- 
ent terms, as it will lay the foundations for our numerical 
investigations that closely follow the loop expansion. 



4.1 Form of the Action 

We are looking for a RG-flow of a given theory defined 
by its generating functional Z . To be definite, let us work 



For introductions into the theory of the renormalization group, see e.g. |34ll4| 
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with the theory of a vector field Vi, and write Z in the 
following way: 



Z = cx-p{-W[{J}]} = J Dvcx-p{-S) 



Di; exp I - - (w„ P^^^ Vj ) - { , Q„^^ vj ) 
-Sint[v,A;Ao] - So[J, A; Ao]} . 



(85) 



(86) 



The action depends on two momentum scales A and Aq. 
By ylo we denote the scale on which we impose the ini- 
tial renormaUzation condition - e.g. the value of the four- 
point-function is fixed to a certain value A4 if all external 
momenta equal Aq. 

The term 5*0 might look uncommon, but is necessary to 
pick up terms nonlinear in J that will be generated by the 
RG-flow. As initial condition, we set So[A — Aq] — 0. 
Starting from a renormalized action on scale ylo, the flow 
is going to generate the renormalized action on all lower 
scaled A, which is the second momentum scale involved. 
From the RG-perspective, S[A = Aq] plays the role of the 
initial condition of the flow. 

The flow equations depend on the choice of the kinetic 
action, so we will define the kinetic term to be 



have been tried. In practice, some propagators will lead to 
simpler rmmerical calculations than others. A very special 
choice of C is the sharp cutoff C^: 

C7'(^)=l for \P\<A, (93) 

C7'(^)=0 for \p\>A, (94) 

which leads to the Wegner- Houghton equation and will be 
treated separately. 

In a similar matter, we define the following kinetic terms 
for anti-commuting Grassmann variables: 



(95) 



where ^ is an antisymmetric matrix in the indices /i 
and v. 

For completeness, we already mention here that we will 
expand the interaction part of the action, 5'int, in powers 
of the fields to illustrate some examples. 



Sint,k, 



(96) 



S'kin = ^iv,P ^w) 



(87) 



where we will call 



1 f ,D ^ . f ^p-l fP^\ . ^ ^int.fc - / Afe (d^p, T;(p,)) S{Y,Pj) (97) 

^^^-^ J d pd qSip + q)v{p)P [^Jviq), J ^ 



where we define 



P- 



P_ 

A^ 



P_ 

A^ 



(89) 



C is the cutoff-function, which has the following proper- 
ties: 



C 



C 



-1 / P 



-1 / P 



A^ 

2 



A^ 



1 for \p\ 



for IpI — > oo. 



G-\\) 



1 



(90) 
(91) 
(92) 



a /c-vertex. The derivation of the flow equations does not 
depend on this expansion; but it is useful in some definite 
calculations. 



4.2 The Wilson Equation 

4.2.1 Integrating out degrees of freedom / Lowering the 
Cutoff 

A simple way to derive the RGE is to calculate the effect 
of a change of the cutoff on an action, keeping in mind 
that both the generating functional and the correlation 
functions may not change 



W -.^ A 



Though it is by no means necessary, one usually assumes 
that C^^ is monotonous, and that it is a smooth approx- 
imation of the step function, thus suppressing degrees of 
freedom on scales bigger than yl, while not effectively al- 
tering those on scales below. The last equation is am- 
biguous, but we define a value for C^^(l), so that the role 
of A becomes definite. We will say that the degrees of free- 
dom that are suppressed are "integrated out", as this part 
of the involved integrals can be interpreted as already be- 
ing performed. Apart from the properties (|90ll92p . we are 
free to define C in any way we like. It follows that not 
even ([55)) is enforced; other definitions of the propagator 



dW 
~dA 



(98) 



{Jvi,Qv,]vj) + 5'int + 5*0^ 



0. 



(99) 
(100) 



In our case, we will lower the cutoff by lowering A, leaving 
Aq as a unit of measurement unchanged. 
Applied to the vector theory, for example, we arrive at the 
following equation for the interaction term of the action: 



5'int — 



1 



Svj 



Sv 



A nice derivation of the Wilson-flow equation is e.g. found in |2]. 
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where we dropped a field-independent term. Taking the 
kinetic term into account, we find the simple equation 



(102) 



The term 



2 Ip 



ss t> 



ss 



will 



the link-term of the flow-equation, while we will call 
~1 1 W'^'Vji §§' loop-term. These names will be justi- 
fied in the following subsection. In complete analogy, equa- 
tions for theories involving Grassmann variables -0* and 
^ with propagator (j95p can be derived: 



This graph gives a contribution to the (rii + n2)-vertex, 
proportional to 

J A„i+iA„2+iP S{Y^p^+p)6{Y^qj -p)d°p. 

(105) 

The two (5-functions will pose problems in the numerical 
process, as they will have to be approximated in order to 
be suitable for a derivative expansion. We can eliminate 
one (5-function by realizing that 

from now on be called Si^^pi +p)6(y^^ qj — p) = Sj^^pi +p)5(^^pi + qj), 



5Vi 



5* 



5S 



5S 
5W„ 



5^u 



5S 



5S 



P.7/ 



(104) 



In case of anti-commuting fields it is important to keep 
track of all extra signs that arise. 

These equations describe the lowering of the cutoff, or in- 
tegrating out of degrees of freedom. Before we proceed, 
let us remark on the graphical interpretation of the RG- 
equations. 



(106) 

which gives us at least the overall conservation of momen- 
tum that we need. 

The loop-term is equally easy to understand: From a ver- 
tex with n + 2 attached lines, two are joined by a prop- 
agator P (Fig. 2]). Again, the symmetry factor is given 
correctly. 




4.2.2 Graphical Representation 

Let us begin with the interpretation of the link-term. For 
the time being, we assume that it is applied to a part of the 
interaction term of the form (^7)) . a vertex with ni -I- 1 at- 
tached lines. Then the functional derivative of this gives us 
a vertex with ni lines; the missing line is linked by the part 
of the propagator that is integrated out, P, to a second, 
similar vertex with, say, n2 + l lines. The graphical result 
is shown in Fig. [31 Observe that the functional derivatives 
automatically lead to the correct symmetry factor of the 
graph. 




Fig. 3: Link-term of the Wilson equation in its graphical 
representation. 



Fig. 4: Loop calculated in the Wilson equation. 
This graph gives a contribution to the (ri)-vertex; 

J A„+2P(^^j<5(^p.)rfV (107) 

The overall (5-function does not pose any problems in this 
case. 

So far, we explained the effect of the flow equation as only 
S'int is concerned. Let us now investigate the contributions 
of the kinetic term. 

The loop-term generated from the kinetic term (Fig. [5]) is 
trivial , as it is field-independent and can be dropped. 




Fig. 5: Field- independent loop that is constructed from 
the kinetic term of the action. 
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Let us consider the terms arising when one field derivative 
in the link-term acts on the interaction, and the other one 
on the kinetic term; this one is compensated by another 
term in the RG-eq.: 



Svi 



'55'kin p SSint 



Sv,, 



(108) 



The remaining term to be considered is 



2./p Sv 



6v 



1 



= — / vP-^PP~^v. 



We can use 



(109) 



(110) 
(111) 
(112) 



and this is, as defined in (|88p . the change in the kinetic 
term. 

Let us summarize: We have seen that the RG-fiow can be 
expressed graphically. Iteratively, we calculate the contri- 
butions from all graphs with one propagator P, and all 
other propagators P, that are inner propagators which 
have been generated in RG-steps before. This is simply 
an application of the product rule: 



i = -A- (PP-^) ^ 
=^ -\ j vP-^PP-^v = ]- I vP-^v, 



(113) 



If we formally define Q[f{p)] to be the sum of all possi- 
ble Feynman-graphs of the theory with inner propagators 
/(p), we can even write down a formal solution to the 



Wilson equation: 

S[A-Aa] = S[A = Ao] 



S[A ^Ao]- I g 

I A 



An 



P(^) 



dA 

-T (114) 
A 



Ap(£) 

dA 



S[A = Ao]-g 



In the limit An ^ oo this becomes 



S[A; Ao] ^S[A^Ao = oo]-g 



dyl(115) 
(116) 

(117) 



This reveals the meaning of the changes to the action: The 
RG-fiow sums up the part of the propagator that is cut 
off iteratively. Notice that we seem to subtract all graphs 
- this is because we are working with e""^ rather than e'^. 
Of course the derivation and graphical interpretation of 
the fiow-equation does not apply directly to the case of a 
sharp cutoff as defined in (|93l IM)) . As our numerical ap- 
proach favours the Wegner- Houghton equation, we discuss 
it in the following paragraph. 



4.3 The Wegner-Houghton Equation 

For numerical purposes, it is easiest to work with the sharp 
cutoff function Cg defined in ([Ml This changes the 
form of the fiow equation drastically. Again, we will not 
present a derivation here as it is found in the original lit- 
erature |31j . but present a short overview. We start from 
dividing degrees of freedom into a high-momentum part 
that is to be integrated out, and a low- momentum part 
that is kept. Expressing the integrated part of the func- 
tional as a change AS to the action, we find 



= J Vvcxp{- J L}, 



(118) 



where the prime denotes integration over the momentum 
shell between A — dA and A. We now expand the action 
two second order in the fields and get 



ss I s^s \\ 

if 5^S ^ SS SSfS^SY^5s\ ISS fS'^SY^ 6S]] , 

where we completed the square. Integrating out the field v in the shell of momentum, which is assumed to be of the 
thickness ^ ^ 1, we get 



e"-^-^ ^ j Dv exp < 




= j Dv exp ' 





AA f , S'^S\ ^ ,f'lSS/'5'^S\ ^ SS 
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and with the aid of 



(det Ay ^ exp{a In det A} = expjaTr In A} (122) P2 



we can finally write down the Wegner-Houghton equation 



S 




5v 5v \ 6v'^ 



- Trln 



(123) 




In the derivation it is used that it is sufficient to work in 
one-loop-order, as higher order contributions are also of 
higher order in dA. A proof of this is found in the original 
work 1311. 



Notice that (|123p seems to differ from the Wilson equa- 
tion p03p by an overall-sign; but this is explained as P is 
negative. 



Eq. (|123p is especially convenient for numerical applica- 
tions, as the contributions to the integrals can be calcu- 
lated explicitly. The logarithm looks problematic at first 
glance, but we shall demonstrate in the next paragraph 
that it has a very simple graphical interpretation, and is 
thus favourable for our graphically based program. 



4.3.1 Graphical Representation 



Once again we start with the interpretation of the link- 
term, which for the Wegner-Houghton equation does 
not involve the bare propagator alone, but the quantity 

(j^S"^ ■ '^liich is more than just the inverse of the ki- 



netic term. Using the geometric series, we can write 



'kin 



kin 



6v^ 



' Jkin 



kin 



' Jkin 



(124) 
(125) 

(126) 



This is the sum of all graphs with i vertices, linked into 
a line by propagators. The first and the last vertex in the 
line are also attached to propagators, which link them to 
terms |^ . Again, let us assume first that these act on the 
interaction part of the action, thus the chain described 
above is linked to other vertices. We therefore find the 
graphical representation Fig. El 



Fig. 6: Link to be calculated using the Wegner-Houghton 
equation. This graph gives a contribution to the vertex 
with Hi outer fields. 



In a similar way, we derive the graphical representation of 
the loop-term. We re-write the logarithm as: 



Trln 



Trln 



Trln 



Trln 



(129) 

The first term is field-independent and is dropped. The 
second logarithm is expanded as a Taylor-series, reading 




(-1) 



4 + 1 



'<^'^k 



6v^ 



(130) 



The corresponding graph is shown in Fig. [T] It is similar 
to the link-term before, but closed to a loop by the trace. 
The factor ^ compensates the rotational symmetry of the 
graph. 




Fig. 7: Loop calculated in the Wegner-Houghton-equation. 
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As before, we still have to sort out the link-terms involving 4.4.1 Field Strength Renormalization 
the kinetic action. Let us start with an example. 



ri r2 




P- 



Ql 



Q2 



Fig. 8: Graph of the Wegner-Houghton flow, linking a ver- 
tex to two outer propagators in this case. 



The graph in Fig. [5] is obviously one of those that arise 
from the link-term; the reader may focus his attention to 
one of the PP~^ legs. Integration is again over the mo- 
mentum shell, so PP~^ = 1. The result looks like the 
vertex, but with the difference that we know that fields qi 
and q2 are depending only on momenta less than A — dA. 
We can thus interpret these terms as integrating out the 
momenta on remaining fields. Integrating out more outer 
fields at the same time would again be of higher order in 
dA, and is thus omitted. 

The change in the kinetic term itself is again simple, and 
not even a sign problem arises as in the Wilson-case. The 
corresponding graph is 



P 



-1 



P 



P 



-1 



We started our integration step with the kinetic term 



'S'kin — 



c^{p)4>{q)KP + <l)- (133) 



We thus defined S'kin to be the only term quadratic in 
fields and momenta in the limit p — >■ 0. After the integra- 
tion step (lowering the cutoff from Aq to A), new terms 
that, according to our definition, should belong to the ki- 
netic term, can be found within the interaction part of the 
action. Such a term has to be counted as kinetic term by 
hand. In practice, the first contribution to the kinetic term 
arises in the second step of the RG-flow as it is of two-loop 
order. The simplest graph contributing to field strength 
renormalization is the so-called sunset graph. Fig. 1101 




Fig. 10: Sunset graph, leading to the simplest contribution 
to the field strength renormalization of 0^-theory. 



Fig. 9: Graphical representation of the change of the prop- 
agator in the Wegner-Houghton-flow. 



and as P ^PP ^ = P ^, this is exactly 

dp vP~^v. 



(131) 



Again, this is precisely the change of the kinetic action, as 
expected from Eq. 

As in the case of the Wilson equation, we are now able to 
give a formal solution to the WHE. The final result reads 



S[A;Ao] = S[A^Ao] 





■ 1 " 


/ ^ 

1 A 


p2 



d^p. 



(132) 



4.4 Kinetic Action - Renormalization and Rescaling 

The field renormalization step is completely artificial. Let 
us clearly discuss the renormalization step, as it seems to 
have lead to confusion. For definiteness, we will demon- 
strate the concept using 0'*-theory in _D-dimensions; for 
other theories the procedure works in exactly the same 
way. Let us emphasize that this step is not unique to the 
ERG, but also found in perturbative renormalization. 



In our case, a graph analogous to Fig. \W\ is to be com- 
puted by our numerical approach in an iterative way later 
on, summing up contributions from every infinitesimal in- 
tegration. The result depends on the used renormalization 
scheme; by means of a Taylor-expansion, we are always 
able to identify the contribution to the kinetic action; we 
will call this 



d^p 
(27r)^ 



xC- 



P_ 

A^ 



d^'q 

/ <P{p)(l}{q)S{p + q)- 



(134) 



This term is transferred from the interaction term to the 
kinetic action by hand, which then becomes 



'-'kin 



1 



xC- 



d^p 

A^ 



q 



D 



(2^) 



p^m<t>{q)5{p + q). (135) 



As our renormalization condition for the field strength, 
we will choose that the coefficient of the kinetic action in 
the limit p — >■ is equal to i. According to the definition 
of the cutoff-properties, we introduce the field-strength 
renormalization factor Z in a way that compensates for 
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the new term in 5'^i,j. If we write for the original action 
((1331) 



Skin — Z- 



d^p \ f d^q 



(27r)^ ; V(27r) 
xC-l (^^^p^cj,{p)q)iq)5ip + q), (136) 

we conclude that Z transforms as 

.l^lnZ=-,y. (137) 

(The sign changes as in the integration step, we actually 
lower A.) The initial condition has to be 

Z[A = Ao] = 1, (138) 

so that (|133p is fulfilled. We can integrate (|137p easily to 

'AoV 



A 



(139) 



Z is now absorbed into the fields in the following way, 
which explains the name of field strength renormalization: 



Ao 
A 



(140) 



We will use this in the next paragraph to determine the 
scaling of the field. This is in complete accordance with 
renormalization conditions met in perturbative renormal- 
ization, see for example [28, 7J. 
The kinetic term now reads 



S' - 



d^p \ f d^q 



xC-i (^^) p2 + q)^ (141) 

which, expressed in renormalized fields, is exactly of the 
same form as (jl33p . 

The point stressed by Golner [16J and Bervillier [6J is that 
we have to do the rescaling step very carefully, to account 
for the renormalization step consistently. 



4.4.2 Rescaling 

The last step in the renormalization group process is the 
rescaling of the momenta, and the functions thereof. We 
define new momenta p in the following way: 



P 



P, 



so we replace 



P^iTolP 



(142) 



(143) 



We can directly see that this replacement also changes the 
cutoff function in the expected way: 



(144) 



so we are ready to identify the new with the old cutoff. 
From the renormalization step, it is now easy to deduce 
the scaling of a field. Beginning with the original kinetic 
term p33p . we conclude that, as S'kin does not scale at all, 
the rescaled field 6 has to be scaled as 



m - ( ^ ) Hp)- 



We call the canonical dimension D^^can 
tells us the anomalous exponent: 



(145) 



Eq. (Hia 



4>'{p) 



A^ 
A 



(146) 



We will now see that due to the effect of the rescaling of 
the cutoff-function, the renormalized kinetic action indeed 

2 

does not scale. As the cutoff is a function of the ratio ^ , 
a change in p has the inverse effect as the same change in 
A. By this, we reverse the effect of the integration step. 
As a part of this, terms are re-distributed back to the in- 
teraction of the theory. Now scaling back the kinetic term: 
The canonical scaling of the fields is compensated by the 
integration measure, and the anomalous scaling by the 
renormalization step. We thus trivially find 



'-'kii 



d^p 



^2 



d^q 



xC*"' (^) P' i'iPWimP + 9), (147) 



again formally identic to p33p as a function of the rescaled 
quantities, as desired. We will drop the tildes and primes, 
formally getting back to (|133p . 



4.5 The Interaction Terms 



As we have explicitly discussed all steps when applied to 
the kinetic action, we now just have to apply them to the 
interaction terms. Let us, as an example, have a look at a 
four- field- interaction 



d^'pi / d^P2 / rf^'ps A4 X 



X(t){pi)(t){p2)(j}{pz)(t){pA)5{pi +P2 +P3 +P4)- 



(148) 



4.5.1 Renormalization 



Once we change (f> to the renormalized field 0', without 
changing the interaction term, we have to renormalize the 
coupling as follows: 
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Pi 



P2 



d^P2 



d^P3 Xicf){pi)(l){p2)(j){pz)(t){-Pl -P2 -Ps) 



d^Ps A4 



^ ) '/''(Pl)0'(P2)'^'(P3)'^'(-Pl -P2 -Ps) 



A4 — A4 



<^ A4 



We need the change in A for an infinitesimal integration 
step, which is 

(149) 



A4 = -4|A4, 



or, for a general interaction Si^t with any number of ver- 
tices, 

Sint, Ren ^ j ^^^^^ 

as the operator j 4>j^ counts the number of fields in a 
vertex. 



4.5.2 Rescaling 

This step is now straightforward. Let us sum up the con- 
tributions: 

— Integral: For each integration measure, we get a fac- 
tor £), and there is one integration measure less than 
there are fields (because of the (5-function) , so we get 
a contribution 

^int,dp j 0^ + ^5int (151) 

(we are still lowering the cutoff). 

— Momentum: The vertex could and will depend explic- 
itly on the momentum, so we introduce another op- 
erator J 4>{p)p s,^{p) I'^^l counts the powers of 
momenta in each vertex. The prime at the derivative 
indicates that it is not acting upon the momentum 
conserving J-function. We get the contribution 



(j}{p)p 



(152) 



^dp J S(j>(p) 

Fields: As derived above, each field brings a contribu- 
tion proportional to — 



Si 



D + 2 



int,0 



so in total we find 

V f ,^'5'int 



(153) 



— Renormalized Coupling: Not to be forgotten, we renor- 
malized the couplings according to (llSOp . so it scales 
itself anomalously, exactly compensating the anoma- 
lous scaling of the fields: 

So, to sum up all contributions, we finally get the rescaling 
term 



int, Rescaling 




(155) 



4.6 The RG-Equation 

For the interaction term, we thus find the fiow equation 



'S'int — 



'^S'int p SSint 




(156) 



This equation clearly depends on the choice of propagator 
([59)1 . To give an example: If we had defined the inverse 
propagator to be 



(157) 



as often found in literature, we would have ended up with 
the same equation as Bervillier [6J or Golner [16J, namely 




(158) 



which in turn is completely equivalent to Wilson's equa- 
tion. 

We want to point out that both equations p56l) and 
(|158p are in a way correct, even though they seem to dif- 
fer by a sign. The equivalence is obscured by a different 
choice of propagator functions, taking advantage of the 
reparametrization-invariance of the equation. 
Propagator (jl57p indeed has some advantages, as in prin- 
ciple other values for fc are also possible, and simplify the 
implementation of K41, as we will see. On the other hand, 
the derivative of (|157p is a bit uglier, and especially the 
derivative expansion gets mixed up. 

We can now aim for the final assembly of the RG-equation, 
written down for vectorial and Grassmannian fields: 
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2 .L ^1 Sv,. 5v,, 



dvi dvi dv.. 



6S SS 



S 5S 



ss_ 

Hi 



2 Jp oi>i 
&_5S_ _ /■ ,* d' 5S 

-DS. 



, d' SS 



(159) 



This equation is general enough for all our applications, 
including different ways of considering the functional de- 
terminant Eq. (j35p . 



5 Derivative expansion 

Clearly, any action has to be approximated to be of nu- 
merical use, as it represents infinitely many degrees of free- 
dom. A common way of approximation is the derivative 
expansion, see e.g. |17| and |26) . In principle, applied to 
the scalar theory, it amounts to expanding the action in 
powers of derivatives: 



Vi^{x))+Oid'), (160) 



in contrast to an expansion in powers of fields, which can 
be seen as expansion around a weak field. As a special case, 
in the Local Potential Approximation (LPA) the action is 
reduced to an interaction term depending only locally on 
the field (f>{x) (and not on its derivatives) and a kinetic 
term whose coefficient Z is held constant throughout the 
flow: 

SLPA^^Jjdc^f + V{4,{x)). (161) 

If we apply the loop- and link-terms to the action ex- 
panded in powers of fields, we are led to rate equations 
for the coefficients. Let us, as an example, apply the Local 
Potential Approximation (LPA) to Eqs. (|105p and (|107p . 
expanded in powers of fields. Starting with Eq. (jl05|) for 
the link-term, we see that the non-trivial part of graphic 
is proportional to 



I A„+i(pi,...,p„,p)A,„+i((ji,...,g™,p) 

xS{pi + P2 + ■■■+Pn +p)S{qi + q2 + ■■■+qm -p)d^p. 



(162) 



As explained above, we transform one of the (5-functions 
into an overall momentum conservation, and use the other 
to evaluate the integral: 



+ni=P i ^n+l{Pl, ■ ■ ■ ,Pn, - ^Pi)>^r,i+l{qi, ■ ■ ■ ,qni,^qj)SC^Pt + ^ 



<l3J 



(163) 



In the LPA, we approximate the couplings to be 
momentum-independent, and develop the cutoff to zeroth 
order in the momenta, so: 

An-i-m,LPA = bm P j—^ A„+i,lpaA„i+i,lpa- (164) 

p^o \ A'^ J 

A difficulty is the momentum-dependence of the factor 
P ^ ^2^' ^ which we need to expand, according to the 
derivative expansion. The result obviously depends on the 
choice of the cutoff; if we apply it to the LPA, we can 



summarize the result into the constant 

Pi := lim P (^^^\ . (165) 

If the cutoff is an approximation of the step function, we 
would expect the limit to converge, and Pi = 0. This 
clearly is not an option, as it would suppress the non- 
trivial character of the RG-flow. 

On the other hand, the loop-equation (|107p leads in the 
LPA to 
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^" = / ^ ( ^2 ] ^n+2bl, • ■ • -p)S{pi + P2 + ■ ■ . + Pn)d'^ P (166) 



A„,LPA = A„+2XPA j P ) d°P (167) 
= A„+2,lpa/?d-i / f f 4 ) ) p'^dp. (168) 



Again, the integral depends on the choice of the cutoff; for 
the LPA we will write 



dp 



p'^dp ^ Po- 



(169) 



Rather than to specify a cutoff, in the LPA it is sufhcient 
to define the constants Pi and Pq. In higher orders of 
the derivative expansion, additional information concern- 
ing the cutoff will be required. 

In the case of a vector theory in three dimensions, the 
situation is not that simple, as products of the type ViVi 
or any contraction with other three-component fields will 
be present. We need to keep track of this to calculate the 
contributions to a renormalization group flow, so we pro- 
pose to expand the terms of the action in powers of fields 
and momenta in the following way: 



qy(Ai,A2,...Ae 
(xi,X2,x^) (ri,r2,...re) 



x,r.q.A 

X {vv)^- {uu)^^ iff)^^ (vu)^^ (vf)^^ {uf)^<^ 

X {vii}*y' {vipy^ {ui^y {uipy^ {fi'*y' {ftpy 0^*^^)" 

X (01)^1(02)^^(03)^^ 

(170) 

From now on, we will work with the coefficients 



qy{Ai,A2,...Ae) 

(2:1,2:2,2:3) (l-l ,I'2,---I'6) 



(171) 



In first order, the terms of the derivative expansion are 
even more complicated, as we also have to keep track of 
terms like piVi{q). 

As the overall number of momenta is fixed for each term, 
and the action itself is a scalar, we get the following pos- 
sible values for the indices of V: 

x^ e No, A, G No, (172) 

r, G{0,1}, gG{0,...,i:»-ri-r3-r5}, 

(173) 

ri+r^+rz^r2+ri + rfi, (174) 

and for Z equivalently. For completeness, we also have to 
keep track of the number of time derivatives Der in a term 



- so the quantities we will work with are the coefficients 

(l)yl^\Der]. (175) 

For the final assembly of our rate equations, we still need 
the scaling dimensions and exponents k for the involved 
fields. The canonical dimension for the velocity field v is 
derived from the energy flow 77', see |13| : 



n', oc 



oc e. 



(176) 



from which we see that \/Jv^V)) oc 1^ . In wavenumber 
space this implies 



4 
3' 



-4- 



(177) 
(178) 



The scaling dimensions of the non-physical fields and con- 
stants are then straightforward: 



(179) 



[p]^-D+l, [^*] = [^] = -| - i, (180) 



[A] = ^-3' ['^^ 
as are the exponents k: 



D 



1 Vi, 



ky = D + fc^, = 2 Vi, 



kf = D- 



ku^-D 



(181) 

(182) 
(183) 
(184) 



We thus finally arrive at the rate equations that we simu- 
lated numerically; here we only present the LPA: 
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[Der] =- 



D 



Vv 



2 



Pv,l 
(Pu. 



^ Pv,0 
Pufi 



-D 



Vf 



P 



/4 



dVdV\ 
df Of J 



+ w 

2 + ?70 



P^. 



dV dV\ 



(x),(A) 



1 

3 

-D 



2 



(x),iA) 

(2Ai + A4 + A + ri + rs) 



i^)XA)\ 

q.ir) J 

i^)XA)\ 
(x).(A)N 



q2y \ (2;),(-4)N 



?7n 

2 



(2^2 +^4 + ^6 + ^3+ r4) 



-i - ^ ) {2A3 +A5 + Ae+r5+ rg) 
2" 



D 1 rfxp , 
-77+0^"^ + »'2 + ?'3 + ''4 + ^5 + rg + 2(7) 



E 



T 



2 



-Ller + D 



Here we defined, for example, 



(x),{A) 



dV dV\ 



(185) 



as the contribution of -§^4^ with the indicated field ex- 
pansion. 

Again, if the determinant pSI) is dealt with in a way that 
implies additional fields, these have to be included into 
the RG-equation in an analogous way. 
Of course, it is in a way a drawback to expand the action 
in powers of fields and momenta. On the other hand, this 
enables us to work with a very simple and fast numerical 
algorithm that we describe in the next section. At this 
stage, we are able to check the scaling of the two- and 
four-point functions {vv) and {vvvv) near the free fixed 
point: 

d , , 2 , , ~ (vvvv) 

dt 6 [vv) 

Correspondingly, we get for the four-point function: 

d 4 ~ (vvvvvv) 

— (vvvv) = —- * (vvvv) -I- 55 * P„ * — ; — ; — (187) 
dt 3 ' [vv) 



We see that in the limit of small couplings, the two- and 
four-point-function scale as predicted by K41: 

d 



dt W = -3*^""^ 



(vv) - (x) 3 , 



— (vvvv) 
dt ^ ' 



'— * (vvvv) 
o 



(188) 
(189) 
(190) 

^ (vvvv) {x)i . (191) 

(Note tire signs: We lower A, which is equivalent to in- 
creasing p, and again decreasing x.) Moreover, Py^Q can 
now be interpreted as a measure for the coupling, defining 
tire meaning of being near the free fixed point. 

6 Numerical Analysis 

6.1 Choice of Renormalization Group Equation 

As indicated before, we tried different choices for the cut- 
off, including the sharp cutoff of the WHE, Eq.([n21). We 
found that this choice is especially suitable for our simu- 
lations, as it allows us to compute the contribution for an 
infinitesimal integration step independently of the coup- 
Hngs involved. 
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We worked with a predictor-corrector-, as well as a Runge- 
Kutta-integration algorithm, both with self-adjusting step 
width. We used two sets of algorithms - one of them in- 
volves explicitly programmed versions of the rate equa- 
tions, while the others worked out the loop- and link- 
graphs automatically, only needing the parameters of the 
physical system. 

Apart from the algorithm of the flow simulation itself, we 
developed a number of tools for the analysis of the result- 
ing data. As the coupling space, in which we are working, 
is very abstract and high-dimensional, it is helpful to make 
first steps by the study of unphysical toy systems, i.e. sim- 
ple and solvable physical systems, and reduced turbulent 
systems (Burgulence) first, to gain the confidence that the 
algorithm is working correctly, and to develop some intu- 
ition for the work with renormalized couplings. 



6.2 Non-Turbulent Systems Analyzed 

We started our investigations by analyzing unphysical 
(toy-) systems with arbitrary constants and dimensional- 
ity of space, to learn more about the detection and features 
of different sorts of fixed points. A main question was how 
structures in coupling space can be recognized, if the di- 
mensionality of the coupling space is high, and whether 
terms of higher order in the field expansion contribute as 
a correction. 

In a second step, we applied the algorithm to standard 
physical systems with known properties, as the scalar and 
the 0(3)-symmetric field theory. The idea behind this was 
on the one hand to ensure that the algorithm works cor- 
rectly and to see how closely we can reproduce analytic 
values for fixed point scalings; and on the other hand to 
approach turbulent hydrodynamics in a stepwise manner, 
interpreting it as a special case of the general 3-vector- 
model. 



iteratively. Following this procedure, the simulated trajec- 
tories approach the ideal trajectories, i.e. the flows directly 
ruiming into our out of the fixed point. 

The shooting method is limited by the numerical accuracy 
of the computer program, and the stepsize adjustment of 
the flow integration, as the algorithm slows down drasti- 
cally when a fixed point is approached. 

In a simulation involving more than two couplings (as is 
usually the case) the projection of the flow onto a two- 
dimensional subspace will in general not look so evident, 
but quite similar if the fixed point is approached closely 
enough. 



-0.003 - 



-0.004 - 



-0.005 - 



-0.006 - 



-0.007 - 




-0.008 I 1 1 1 1 1 1 1 ' J 

0.11 0.12 0.13 0.14 0.15 0.16 0.17 0.18 0.19 0.2 

Fig. 11: Fixed point of an unphysical model system, as 
found by use of the shooting method. The attractive di- 
rection goes to the upper right and lower left corners of 
the diagram, the other two directions are repulsive. 



6.2.1 Toy Systems 

We worked with a number of unphysical systems for trials 
of the algorithm and analysis tools, thus merely looking 
for interesting structures. These systems were defined by 
an action consisting of a propagator, a two-field- and a 
four-field-interaction, where the field was a 3-vector-field. 
Parameters were deliberately adjusted to allow the pres- 
ence of different fixed points. 

Investigations of the coupling space were mainly done us- 
ing the shooting method, which is especially useful for 
finding fixed points. In practice, one initiates a number of 
RG-flows, starting from initial conditions sufficiently close 
to each other, and searches the topology of the flow for in- 
teresting structures. An example can be seen in Fig. [Til 
where the location of the fixed point is easy to guess. 
To identify the location more precisely, one would re- 
peat the method with initial conditions closer to the 
guessed fixed point couplings, leading to a picture simi- 
lar to Fig. 1121 In this way, one approaches the fixed point 



-0.003 



-0.0035 - 



-0.004 - 



-0.0045 - 



-0.005 - 



-0.0055 - 



-0.006 




0.14 0.145 0.15 0.155 0.16 0.165 0.17 0.175 

Fig. 12: Fixed point as in Fig. [TT] after adjusting the ini- 
tial conditions. Position and ideal trajectories of the fixed 
point can clearly be identified. 
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Other interesting structures include focus fixed points; an 
example can be seen in Fig. [131 As explained by e.g. Sor- 
nette [30], these fixed points can be associated with com- 
plex eigenvalues, and can thus be identified with some 
discrete scale invariance. Observe that the simulated tra- 
jectory ends at the fixed point due to numerical errors, 
rather than to encircle it and approaching it asymptoti- 
cally. 

0.0295 I , 1 1 1 , , 1 




0.0255 I ' 1 1 1 ' ' 1 

-0.135 -0.13 -0.125 -0.12 -0.115 -0.11 -0.105 -0.1 

Fig. 13: Focus fixed point of an unphysical model system, 
as found by use of the shooting method. The repulsive and 
the attractive directions can clearly be distinguished. 



The advantage of our approach is the fast integration of a 
large number of couplings, and in that way circumventing 
the drawbacks of the expansions. Simulations were done 
with up to 100 couplings, though it has to be said that 
the identification of fixed points becomes nearly impossi- 
ble in these high-dimensional spaces. Working with such 
a number of terms can only be done iteratively, meaning 
that one starts with a low number of couplings, identi- 
fies the fixed point and than changes to more and more 
terms, hoping that these act as corrections to the overall 
behaviour. 



6.3.1 Fixed Point Analysis 



pr 
-0.02 - 
-0.04- 
-0.06 - 
-0.08 - 
-0.1 - 
-0.12 - 



6.2.2 Simple Physical Systems 

The renormalization group flows of the scalar theory and 
the 0(3)-symmetric theory in D dimensions have been an- 
alyzed by P. Diiben [5]. By reproducing known values of 
these theories like fixed point locations and scaling (also in 
the e-expansion), we tried to get closer to the much more 
divert general three-vector-theory, and also to check the 
correctness of our algorithms. We found that we are able 
to accurately reproduce the values known from literature, 
to a given order of the e-expansion. Our results will be 
published in a forthcoming article. 

6.3 Local Potential Approximation of Hydrodynamics 

We now return to the analysis of the action for hydrody- 
namics derived above in the LPA. The system is speci- 
fied by the dimensionality of space and symmetries of the 
fields; action merely serves as the initial condition of 
the flow. 

The simulations were performed using two distinct algo- 
rithms: The first one iterating the rate equations derived 
in the previous sections and doing the book-keeping of 
the terms involved explicitly, the second one finding the 
graphs to be computed automatically. The second formu- 
lation turned out to be not only more elegant, but a great 
deal faster than the cumbersome implementation of the 
book-keeping. 



-0.14 1 ' < < ' < < < < ' 1 

1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2 

Fig. 14: Fixed point, as found by use of the shooting 
method. The plot shows the fixed point value of the (vv)- 
coupling, depending on the anomalous exponent. 



0.03 I , 1 , 1 , , 1 , 1 1 

0.025 - 

0.02 - 
0.015- 

0.01 - 
0.005 - 

qI '-' ' ' ' ' ' ' ' ' ' 

1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2 

Fig. 15: Fixed point, as found by use of the shooting 
method. The plot shows the fixed point value of the 
(wTO)-coupling, depending on the anomalous exponent. 
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It is a simple calculation to show that for values rj > 1.5 
of the anomalous exponent, a non-trivial fixed point ex- 
ists in the vicinity of the trivial one. We used the shooting 
method to determine the position of this non-trivial fixed 
point, depending on the anomalous exponent, as can be 
seen in Fig. [T3] and Fig. [TSl The distance to the origin of 
coupling space can be seen to grow linearly with 77, but, 
alas, no physical property could be ascribed to it. Note 
that for 77 < 1.5, this fixed point does not exist. As we 
have already explained, in flow-simulations it is generally 
preferable to calculate j], rather than to search for it by 
means of the shooting method. In the LPA, on the other 
hand, one has ry = as no field renormalization is per- 
formed. 



6.4 Scaling of the Trivial Fixed Point 

It is straightforward to analyze the scaling of the triv- 
ial fixed point. The correlation functions of even orders 
are directly computed by the RG-flow; after Fourier- 
transformation to physical space we can read off the scal- 
ing, and find: 



Order of the Correlation 
Function 


Scaling Exponent 


2 


0,666 ±0,017 


4 


1,338 ±0,035 


6 


1,999 ±0,052 



Table 1. Scaling exponents at the trivial fixed point 



The correlation functions of odd orders are not explicit 
terms of the action and so have to be measured indirectly. 
The correlation function of order n can, for example, be 
derived from the term (uw"), if the scaling of the field u 
is known. We suggest to measure the scaling of u from 
the two-point- function (uu), and subtract it from (uw"). 
In this way the following exponents can be measured: 



Order of the Correlation 
Function 


Scaling Exponent 


1 


0,3334 ±0,0018 


3 


1, 0004 ± 0, 0012 


5 


1,6681 ±0,0012 


7 


2, 3348 ±0,0012 



Table 2. Scaling exponents at the trivial fixed point 



As these numbers reveal, we can conclude that the trivial 
fixed point represents the scaling of the K41-theory. 

7 Conclusion / Outlook 

We have shown how to define a generating functional for 
hydrodynamic turbulence, including a strict treatment of 



the incompressibility condition. The non-local interactions 
have been transformed into local ones, by a re-definition 
of unphysical fields and constants, and by introduction of 
new fields. We have shown how to approximate the result- 
ing action within a derivative expansion. 
On the side of the renormalization group, we discussed the 
procedure of renormalization and rescaling in detail in or- 
der to clarify some obscure point in the literature. Our 
work leads to a RG-equation for the turbulent action, and 
a set of rate equations after the application of the deriva- 
tive expansion. These rate equations have been simulated 
numerically, and the resulting data have been analyzed. 
As will be published in another article, we tested the nu- 
merical algorithm by reproducing known values for non- 
trivial scalings of the scalar theory in 4 — e dimensions, 
and the 0(3)-symmetric field theory. The results are in 
agreement with values found in literature, so the numeri- 
cal algorithm can be considered to be reliable. 
Our numerical work involves the computation of the RG- 
flow for products of Grassmannian variables, as well as 
3-vectors that have to be kept track of according to their 
combinations. 

We were able to identify the trivial flxed-point with the 
scaling exponents predicted by the K41-theory. 
On the other hand, we have to admit that we were not yet 
able to reproduce intermittent exponents for the structure 
functions of turbulence that would agree with the exper- 
imental values. The reason clearly is the complexity of 
the general 3- vector-model, including all theories that are 
based on hydrodynamics. Although the basic foundations 
of these theories are easily understood, all of them (in- 
cluding Navier-Stokes and Burgers turbulence) involve the 
same dimensionality of space and symmetry of the fields, 
while leading to different predictions for the intermittent 
exponents. 

Moreover, it is questionable whether the analysis of a fixed 
point will eventually lead to an understanding of inter- 
mittency. Data seem to hint to a cutoff-phenomenon; the 
probability distribution of the velocity increment looks, 
for small distances, like a Levy-distribution; on large scales 
like normally distributed [H]. We are strongly convinced 
that this can be explained by a crossover between two 
fixed points, but cannot yet justify it by our simulations. 
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